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1. Introduction 

One of the models of accretion discs in a close binary system is the spiral 
shock model. This model was first proposed by one of the present authors 
(Sawada, Matsuda & Hachisu 1986a, b, Sawada et al. 1987). A number 
of authors have confirmed since then that spiral shocks appear in two- 
dimensional discs. In the case of 3D, Sawada & Matsuda (1992) used a 
Roe scheme to calculate the case of 7 = 1.2 with mass ratio of unity. 
They found the existence of spiral shocks, but their calculation was done 
up to only half a revolution period. Yukawa, Boffin & Matsuda (1997), 
who made numerical simulations using SPH with mass ratio of unity, also 
demonstrated that spiral shocks existed in the case of 7 = 1.2. 

The purpose of the present paper is to perform 2D and 3D finite volume 
numerical calculations with higher resolution up to enough time to confirm 
the existence of spiral shocks. 

2. Calculations and Conclusions 

We considered a binary system composed of a mass-accreting primary star 
and a mass- losing secondary star. The mass ratio of the mass- losing star to 
the mass-accreting star was one. We set a numerical grid centered at the 
primary star and just touching the inner Lagrangian point LI. In addition, 
we used Cartesian coordinates having 200 x 200 grid points in 2D and 200 x 
200 x 50 grid points in 3D, respectively. Finally, we assumed a polytropic 
gas with constant specific heat ratio 7. 

In this paper, we examined the cases of 7 =1.2, 1.1, 1.05, and 1.01. 
Symmetry about the orbital plane was assumed in our 3D calculations. 
We used the Simplified Flux vector Splitting (SFS) scheme (Jyounouti et 
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al. 1993; Shima & Jyounouti 1994). In 2D calculations, we found that the 
smaller 7 was, the more tightly the spiral wound. The density contours in 
the orbital plane in the cases of 7 = 1.2 and 7 = 1.01 of our 3D simulations 
are presented in Fig.l. These figures show the existence of spiral shocks 
in the cases of 7 = 1.2 as well as in all the other cases. However, they 
do not show as clear a difference in spiral structure as they do in our 2D 
calculations. (Other results can be seen in CD-ROM). 
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Figure 1. Density contours in the orbital plane at t = 27. Left: the case of 7 = 1.2, 
Right: the case of 7 = 1.01. 

Steeghs, Harlaftis & Home (1997) found the first convincing evidence of 
spiral structure in accretion discs observationally. We have now succeeded 
in reproducing their results using our numerical calculations. (Please see 
Matsuda et al. in this volume). 
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